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Abstract 

We demonstrate that non-exponential decays of unstable systems can be 
understood as yet another example of nonextensivity encountered in many 
physical systems and as such can be characterized by the nonextensivity pa- 
rameter q. 
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It is known that for the unstable system, for which the probability 7 to decay in 
a unit of time is constant, its survival time probability is given by simple exponential 
law 

P(t) = exp(- 7 t). (1) 

This means that our system is not really aging but with passing time shows its 
persistency to decay. On the other hand, it is also widely known [0 that spontaneous 
decays of quantum-mechanical unstable systems cannot be described by the pure 
exponential law (neither for short nor for long times), i.e., one observes that for 
them P(t) > exp(— jt) and, in practice, P(t) exhibits the power-like behaviour of 
the type §Q 

P{t) oc r s . (2) 

We would like to concentrate in this letter on the observation that there is enor- 
mous variety of different physical phenomena which do exhibit a power-like be- 
haviour instead of naively expected exponential one. They can be described in a 
most economical and adequate way in terms of Tsallis nonextensive statistics [[| 
by introducing only one new parameter q, known as the nonextensivity parameter. 
The common feature of these phenomena is that they are characterized either by 
long-range interactions and/or long-range microscopic memories or by space-time 
(and/or phase-space) (multi)fractal structure of the underlying physical process H, 
or by intrinsic fluctuations of some parameters present in the system under consid- 
eration HQ. In this approach one parameter power-like formula replaces the usual 
exponential (Boltzmann-Gibbs) one[|: 



1 - (1 - q) 1 

T 



L q=l {t) = - • exp 

r 



(3) 



The question we would like to address is whether non-exponentiality of decays men- 
tioned above can be expressed in the same approach, by deviation of the suitably 
defined parametr q from unity. 



In we have demonstrated that parameter q (and power-like distribution it 
characterizes) occurs whenever parameter 1/r of the original exponential distribu- 
tion fluctuates according to gamma distribution. In this case q is given by the 

lr The ability to observe such nonexponential decays j^j and interest in quantum version of the 
Zeno paradox || resulted in the increase of interest in the problem of unstable quantum systems. 

2 This last point, which we shall explore here, is receiving growing attention recently |Q. 

3 We would like to stress at this point the fact that eq. (||) is the exact solution of a simple 
nonlinear ordinary differential equation of the type: ^| = ^jr- with y(t) — jjgjj as discussed in 
detail in cf. 0. 
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relative variance of such fluctuations. We shall demonstrate now, using random ma- 
trix approach || that, indeed, the non-exponential decays can emerge in a natural 
way from the possible fluctuations of parameter 7 = 1/r in exponential distribution 
dH). To this end let us start with discussion of some consequences of quantum chaos, 
which emerge from applying random matrix hamiltonians, elements of which follow 
some random probability distribution ||. Namely, let us consider situation where we 
have N intercoupled states with couplings between them described by hamiltonian 
given in terms of N x N hermitian (and rotatiotionally and time reversal invariant) 
matrix (with real matrix elements): H = [Hij] NxN , (Hij = Hji). Suppose that one 
is interested in ensemble of such random matrices with P{H) being a probability 
of a given matrix H. There is a number of propositions how to define P(H) corre- 
spoding to different models of randomness of matrix elements used ||. However, in 
the limit of large matrices (large N) which are hermitian and real, all of them lead 
to results analogous to those obtained by assuming gaussian ensemble^. This can 
be seen as follows. Let P(H) be such that: (i) - elements are uncorrelated, i.e., 

N N 

p{H) = nnw. w 

t=ij>i 

(ii) - distribution (£|) is invariant in respect to arbitrary unitary transformation of 
the basic states (because initial basic states \i > used to define have been chosen 
in arbitrary way). Let us now look for distribution P(Ci) of the amplitudes Cj for 
eigenvectors of H assuming that our hermitean ensemble is invariant in respect 
to the orthogonal change of the basic vectors. Such ensemble is distributed in the 
probability space with some density P(C\, . . . , Cn) with Cf + C% + C% + - ■ -+Cj^ = 1 
and joint distribution function for N amplitudes has the form 



n 



N 



p(c u ...,c N ) = ^-8ii- yen (.-)) 



(where normalization constant Qn is equal to the volume of N- dimensional sphere 
of unit radius and factor 2 originates from the fact that both positive and negative 
Ci enter P{C\, . . . , Cn))- Integrating @ over all but one coordinates Cj gives the 
following distribution P{Cj) of amplitude Cj 



JV-3 



4 The gaussian ensembles were first discussed by E.P.Wigner [|T^|. For analysis of the gaussian 
ensemble considered here sec Ol. 
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which in the limit of large N has the following form [] 



where 



C 



N 
— exp 



1 

N 



N r2 
2* 



and 



exp 



c? 



(Cff) - (Cff 



2 <C?>, 



(7) 



However, for our purpose we are really interested in decay width ji, which is pro- 
portional to the probability Cf of finding single component \i > in the composite 
state. Its distribution is given by 



P( 7i )d 7l = 2P{C l )\dC i 



(9) 



where factor 2 reflects the fact that both positive and negative values of Ci contribute 



to the same value of 7$. Because \dCi 
such widths: 



d-, 



dji one gets following distribution of 



P(ji)dji = 2P(Q 



dji 



exp 



dji (10) 



V27T7i < jt > ~" r \ 2< Ti >_ 
z) = (1 — i2 < 7i > z)~ x l 2 . In the case where 



d%/ dQ 

characteristic function of which is 
contribution to the width 7 comes from many channels, distribution we are looking 
for is obtained by summing all contributions from the particular channels of widths 
7j. For v different independent channels with the same mean values 70 =< 7i > 
the corresponding characteristic function will be (f> u (z) = (1 — i2^Q z)~ u ^ 2 (with 
< 7 >= z/70), i.e., the corresponding distribution function will be gamma function 
of the form 

ft(7)*T = (^—) f^r' exp f-^M • (ID 



r(j) V 2 <7>/ V 2< 7>, 



2 < 7 > 



(12) 



This completes our proof that, as a results of composition of contributions from the 
v channels, we obtain distribution of the widths of the decays ^(7) in the form of 
gamma distribution, fluctuations of which are characterized by the relative variance 
of the form 

((7- <7>) 2 ) = 2 
< 7 > 2 v 

and they decrease with growing v Q 

5 By using Stirling formula for (N >> 1) leading to lnP(Ci 
noticing that for large N we have Cf << 1. 

6 The (O) is for the integer values of v encoutered here identical to the \ 2 distribution obtained 



±ln^ 

2 111 2tt 



/V 



2 ln(l-Cf), and 



in §. 
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We have therefore demonstrated that, indeed (at least in the above example), 
the fluctuations of the life times r = I/7 around mean value Tq = 1/ < 7 > are 
described by gamma distribution of the form ( [TT]) of the variable 1/r. According to 
our philosophy such fluctuations of parameter 7 in the exponential distribution ([!]) 
lead therefore directly to the Levy distribution of the form || [?J 



Tj \T / V TJ Tq 

1 2 



To 



l-q 



(13) 



with nonextensivity parameter q equal to q 



We shall close with some remarks on chaotic systems and mixing of configu- 
rations. Notice first that distribution P v {pf)d^ is known for v = 1 in the nuclear 
physics as the so called Porter-Thomas distribution fl2| f\ and describes in a satis- 
factory way the widths of resonances observed experimentally |14[]. On the other 



hand the many channel case of v > 2 corresponds, for example, to the case of total 
widths for gamma emision or for overbarrier fission. Model using random matrices 
corresponds to the extreme case of total mixing between allowed degrees of free- 
dom. Such mixing of configurations is realized in the case of compound nucleus 
(where energy of excitation is divided among many degrees of freedom). It means 
that component representing one-particle motion in reference to the whole nucleus 
is decomposed into wave function of a large number of different resonance states. 
In the Fermi gas model the compound nucleus corresponds to a particularly strong 
mixing of different configurations, both corresponding to the stationary and to the 
resonance states. Such emergence of the mixing of configurations, which is leading 
to the compound nucleus, can be understood as effect of an increase of importance 
(paralel to the increase of the excitation energy of the system) of small perturbations 
in the one-particle motion. It happens because in this case the number of exciting 
particles increases (as increases the number of possible interactions) and excitation 
energy per single particle increases as well (the density of final states accessible in 
the collisions increases) |. 

One should stress also that random matrix approach is applied in the chaotic 



quantum systems [15] and it seems to be the characteristic feature of such quantum 



7 Deviations from this distribution, which can be attributed to small N involved, have been 
discussed recently in |l3| . 

8 In this respect there is a whole list of problems which could possibly be considered, for example, 
the problem of mixed states in the quantum physics, connection with the composite biochemical 
systems for which fluorescence time distributions are known to be non-exponential (and are usually 
described by a composition of a number of exponents). 
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systems, which in the classical limit show chaotic behaviour. It is connected with 
the fact that in such systems there is no symmetry, i.e., there is no degeneracy and 
no selection rules which would exclude interactions of some particular states f\. 

We are grateful for fruitful discussions with Prof. C. Tsallis. The partial sup- 
port of Polish Committee for Scientific Research (grants 2P03B Oil 18 and 621/E- 
78/SPUB/CERN/P-03/DZ4/99) is acknowledged. 



9 For application of random matrices to description of chaos in quantum systems see [pi 
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